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Abstract. The paper reports on a study of electron-phonon interaction within a 
limited nanosized region. We invoked the modified Frohlich's Hamiltonian to calculate 
the electron self-energy, as well as the elastic and inelastic scattering cross sections. 
New effects have been revealed, more specifically: a bound state forms within the 
limited nanosized region, electrons undergo resonant elastic scattering, with strong 
inelastic scattering being possible from this state even at low electron energies. The 
effect of scattering on the magnetic-ficld-indepcndent dephasing time, in particular, 
in a diamond-decorated carbon nanotube, has been determined. The effect of strong 
inelastic electron scattering on thermal resistance at the metal-insulator interface is 
discussed. 
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1. Introduction 

Interaction of electrons with vibrational modes, phonons, in nanosized regions has been 
attracting recently considerable interest pQ . This interaction underlies many mesoscopic 
phenomena; to cite an example, it was found to affect strongly the current flowing 
through a nanosized region [2], and it is this interaction that mediates the electron 
dephasing time [3j. It is well known that in the low-temperature domain the Kondo effect 
contributes significantly to the magnitude of the dephasing time. This time depends 
substantially on magnetic field [I]. It was predicted recently [5] and demonstrated 
subsequently experimentally that the vibrational mode affects the dephasing time of 
electrons propagating in an amorphous film [5] or a graphene sheet [7J- Interestingly, this 
dephasing time can be quite large, and it does not depend on magnetic field. Electron- 
phonon interaction mediates transport phenomena in nanostructures as well. We are 
witnessing presently intense research interest in composite materials based on the metal- 
diamond pair [8]. Thermal conductivity of such composites depends substantially on 
the possibility of energy transfer between the electron gas in the metal and the diamond 
phonons. This process occurring in the vicinity of the metal-diamond boundary, it 
appears essential that the theory of this phenomenon takes electron-phonon interaction 
in the nanosized region into proper account. 

All studies in the theory of electron-phonon interaction consider either spatially- 
uniform interaction [H] or interaction at a point [5]. In both cases, the approach 
used invokes either the Frohlich or the Holstein Hamiltonian. It has recently been 
demonstrated that inclusion of the spatial distribution of local vibrational modes affects 
noticeably the properties of the system under study [TO]. The author is, however, 
unaware of any publication which would consider the dependence of the properties of a 
system on the size of the region supporting electron-phonon interaction. This suggests 
the need of tackling the fundamental problem of extending the results amassed on point 
and spatially uniform electron-phonon interaction to cover the case of interaction within 
a local region of finite dimensions. 

The goal of the present study is development of a theory of electron transport 
through a nanosized region with local electron-phonon interaction. The study is 
conducted for the particular case of the impact of electron-phonon interaction in a 
local region on current flow along a carbon nanotube decorated by nanodiamonds, a 
novel material obtained recently [UJ. It is maintained [12] that such interaction affects 
markedly the electronic properties of this material. We are also going to consider the 
effect of strong inelastic scattering on thermal resistance at the metal- insulator boundary 
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2. Model 

2.1. Hamiltonian and the electron self-energy 

We start by constructing a Hamiltonian for systems with local electron-phonon 
interaction. As the starting point will serve the Frohlich Hamiltonian into which we 
shall introduce a coordinate-dependent function describing electron-phonon interaction. 
The Hamiltonian of the system can be written in the form H = Hq + H ep . As usual, 

H = Y^ ZpCpCp + ^q^q a q (!) 

v q 

- is the Hamiltonian of noninteracting electrons and phonons, and 

H ep = J drp(r)V(r) (2) 

- is the electron-phonon interaction Hamiltonian. The Hamiltonians are presented in 
their standard notation. We used here operators of creation and destruction of electrons, 
Cp,c p , and phonons, ai~,a q . The energies of the electron, e p , and of the phonon, hu q 
, depend accordingly, on the momentum p and the wave vector q. We shall assume in 
what follows that Planck's constant H — 1. In the H ep Hamiltonian, p(r) is the electron 
density operator, and V(r), operator of the electron-phonon interaction potential. 

By local electron-phonon interaction we are going to understand here the interaction 
which in a certain specified region is stronger than in the surrounding regions. We will 
assume that electron-phonon interaction exists in such a region only. To discriminate 
this region and identify properly its differences from the other parts of the sample, 
introduce the localization function A(r). This function is of order unity within the region 
of interaction, and small outside it. This may be, for instance, a Gaussian function. We 
shall write the interaction potential in the form of a product X(r)V(r) rather than in its 
conventional notation V(r). This form of the potential underlies all of our subsequent 
consideration. We further write Hamiltonian (J2J) in terms of the Fourier transforms of 
the factors it contains. Integration over the coordinate r yields: 

H ep = J2Hk)Y.P(<l + k ) V (<l) (3) 
fc q 

We next express the quantity p(q + k)V(q) in its standard form through the operators 
of creation and destruction of electrons and phonons [13]. Thus we have 



H ep = J2Kk) 



(4) 



Mq C p+q+k C p( a q + a -q 

k lp,q 

This expression contains the coefficient of interaction of electrons with phonons M q . 
We readily see that for the localization function A(r) = 1 Hamiltonian (J3J) transforms 
throughout the sample to the Hamiltonian of spatially uniform electron-phonon 
interaction [§], while for A(r) = S(r) it becomes the point-interaction Hamiltonian 
[5]. We shall formulate now the electron scattering problem within a region with local 
electron-phonon interaction. We shall consider Hamiltonian H ep as a perturbation and 
perform the calculations as this is done in Ref. [13J. Green's function G(p, e) written 
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in the momentum representation for the complete Hamiltonian H is related to Green's 
function G(p, e) corresponding to the unperturbed Hamiltonian H . This relation can 
be presented in the form G = Go + Go^Go , where E(p, e) - is the electron self-energy 
(SE). The procedure of SE calculation is known to consist essentially in calculating 
the interaction Hamiltonian (jlj) averaged over the electron and phonon states. An 
important note is here in order; in view of our subsequent calculations of the cross 
sections of scattering from a region with local electron-phonon interaction, one can 
conveniently consider, from a formal standpoint, a system combining many regions with 
local-phonon interaction whose concentration is n = 1. In this case, one will have to 
perform averaging over all positions of such regions, an approach adopted in [13J. To 
do this, we transform Hamiltonian H ep to operators in the interaction representation. 
The localization function X(r) being time-independent, the interaction Hamiltonian 
contains the dependence on time only in the bracketed factor in Eq. (jlj). Therefore, 
the procedure of SE calculation for local electron-phonon interaction is similar to that 
for electron-phonon interaction in the spatially uniform case. In this Section, we shall 
restrict ourselves to calculating SE in the first approximation. The exact expression for 
the SE will be developed in the next Section. Thus we have: 



where, S u (p, e) is the self-energy of the electron for the case of spatially uniform electron- 
phonon interaction. In order to calculate the SE of electron-phonon interaction in a local 
region with the help of Eq. (jSj) , one will have to know the form this expression assumes 
in the spatially uniform case, E u (p, e) , which, in its turn, depends on the system being 
considered. 

Indeed, in the particular case of Einstein phonons of frequency cu®:, we have 
9E u (p, e) ~ Q(e — w ), where 0(e) is the Heaviside step function. We immediately see 
that the imaginary part of SE depends on energy while being independent of momentum. 
Substituting this expression in Eq. (jSj), we see that because integration is performed 
over the momentum k, the imaginary part of the SE in the case of Einstein phonons will 
be likewise a stepped function of energy. Analytic continuation of this relation yields, 
just as demonstrated in [T3], a log dependence of the real part on energy. In other 
words, consideration of Einstein phonons in a local region does not introduce anything 
new of particular significance. New effects can be obtained by analyzing the behavior of 
phonons with a linear dispersion law. This is exactly the case we are going to consider 
now. 

2.2. Effective potential in local electron-phonon interaction 

For free electrons, the real part 3?S u (p, e) is constant and does not depend on the 
electron momentum and energy. In other cases, the real part 3ftE u (p, e) does depend 
on these quantities. To isolate qualitative effects, we shall assume that in the case 
under consideration the real part e) is constant too, and that it is equal to So- 

Substituting this quantity in Eq. (j5j), we see immediately that e) is likewise a 




(5) 
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constant. We are going to use this result in the subsequent Sections. However to better 
understand the physical aspects underlying this phenomenon, it will be appropriate to 
recover here Hamiltonian H and subject it to Bogolyubov's transformation. To do this, 
we rewrite the initial Hamiltonian by changing in (j3J) the dummy summation variable 
k to k — p — q, and, subsequently, by transforming operators a q to the new set 

s - Kq)— Y. c i c k 

U Q k 

The sum in this expression, as is well known, is the number of electrons. Since the carbon 
nanostructures such as nanotubes and graphene can be a semimetal or a semiconductor 
[TB] the low temperature free electron concentration is small. Thus, it is possible 
use single-electron approximation. On performing the transformation, we see that the 
Hamiltonian (j3J) of the system can be presented in the form 

M 2 

H = Y,^pCpC p - 2^ A(q)A(fc - q)— -c£ +p c p + 

P kpq U 1 

M 2 

+ E^< a <7 + E A (tf) 2 — + E X (k - P - q)M q c+c p (a q + a± q ) 
q q u i k+p,q 

We shall neglect in what follows the term before the last. This is the shift of 

the origin of energy. The first two terms will be considered as the "unperturbed" 

Hamiltonian for an electron in the region of local electron-phonon interaction. This 

Hamiltonian can be recast in the form of the Hamiltonian for the electron residing in 

the effective potential U(r). In the momentum representation this effective potential 

assumes the form: 

M 2 

U(k) = -2j2Hq)Hk-q)^- (6) 

q U q 

Rewriting now Eq. in the coordinate representation, we can say that the electron 
resides in a potential well with a characteristic depth So = ^ q M 2 /oo q . A potential 
similar to the effective potential U(r) is used in the theory of electron scattering from 
an impurity. The specific form of potential ([6]) is determined naturally by the structure of 
the material under study, which is set by the localization function X(r) and the dispersion 
dependences of the frequency u q and of the electron-phonon interaction constant M q on 
the phonon wave vector q. Nevertheless, one can draw certain qualitative conclusions 
of a general nature on the behavior of an electron in a potential generated by local 
electron-phonon interaction. 

We assume the localization function to have the form A(r) = exp(—r 2 /a 2 ). Our 
subsequent consideration depends substantially on the size of the system under study. 
Since the theory developed below is to be applied to carbon nanotubes decorated by 
diamonds, we are going to study quantitatively here the one-dimensional case. The 
three-dimensional case will be treated only qualitatively. It is known that for an 
electron in a potential well there always exists a bound state, both in the one- and 
two-dimensional case. We are going to denote by n = J2m\eo\/h 2 the wave vector 
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corresponding to the bound state with energy e - F° r the characteristic energy scales of 
the energy S corresponding to the polaron binding energy, and the interaction region 
size a of interest to us here, we obtain a shallow potential well with the energy of the 
level Eq = ma 2 T, 2 l /2h 2 . Equation flB]) suggests that the effective potential U(r,a) tends 
exponentially to zero for both small and large sizes a. Whence it follows that one has to 
take into account the decrease of the well depth itself. We assume that the wave vector 
k scales with the size a as k = maHo/h 2 exp(— a/a c ) . Here the parameter a c is the size 
at which k reaches its maximum value, and the level lies closest to the well bottom. For 
qualitative studies of the three-dimensional refinement is in order here, namely, 

that the bound state in 3D exists only if the characteristic size of the well is larger than 



We see that, while proceeding along the lines of perturbation theory, we used the 
free-electron model (jlj) as the starting point, we came already in the first approximation 
to the conclusion that the system should support under certain conditions formation of 
a bound electron state in a region with local electron-phonon interaction. As we are 
going to see in Sectional the existence of such a bound state will produce a considerable 
impact on the amplitude of electron scattering from a region with local electron-phonon 
interaction. In any case, this "phonon" bound state resembles in many respects the 
bound state in an impurity. The imaginary parts of the SE in systems with impurities 
and with local interaction are essentially different. 

2.3. Imaginary part of the self-energy 

Consider separately calculation of the imaginary part of the SE, which is 

intimately connected with the possibility of emission of a phonon by an electron. To do 
this, we have first to specify the form of the imaginary part of the SE, Q ; E u (p, e), for an 
electron in the spatially uniform case. By the theory of unitarity, it can be related to 
the imaginary parts of Green's functions of the electron, QGq, and the phonon, QD 



We assume here for the sake of simplicity that the electron-phonon coupling coefficient 
M q does not depend on q, M = M q . Although for the acoustic phonons being considered 
here, M q does depend linearly on q, this standard simplification [IB] can be done by 
invoking the first mean value theorem for the integral (J7J). 

For spatially uniform electron-phonon interaction, the calculation is performed 
assuming k = p and e = p 2 /2m. For local interaction, these conditions are not met. We 
have therefore to calculate QT, u (p, e) without making use of the relation between e and 
k. Introduce into this equation the expressions for Green's function for free electrons 
and Green's function for free acoustic phonons having a sound velocity c. We finally 
come to 



its critical value J h 2 /2mTiQ. 







(8) 



L 
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where L - is the region of integration over the phonon wave vector q, satisfying the 
inequality: 

M <£ _ C9< (^i±i)! (9) 

2m 2m 
In the case of spatially uniform electron-phonon interaction, condition Q transforms 
to the well-known condition that an electron can emit a phonon only for p > mc. This 
threshold condition specifies the possibility of phonon emission by an electron. For a 
spatially nonuniform electron-phonon interaction, this inequality has always a solution 
for k < \j2me = p. Examining Eq. ([SJ), we see that when inequality ([9]) is satisfied, 
the imaginary part of SE (j5J) becomes nonzero. Thus, in the local electron-phonon 
interaction region there will always be energy losses involved in the emission of acoustic 
phonons by an electron. This is the main difference of interaction in a local region from 
the case of spatially uniform electron-phonon interaction. To conduct a qualitative 
analysis of the 3D case, we assume that for k < p integration in Eq. (jSJ) runs over all 
phonon wave vectors q, from zero to the Debye wave vector k^. Substitute Eq. (jSJ) into 
(JSJ). Integration performed in the pa — > approximation yields 

3£(e,p) = -£oM 2 . (10) 

Thus, in the case of local electron-phonon interaction, the imaginary part of the SE will 
be proportional to the electron energy and does not have a threshold. The dependence 
of the SE, Eq. (jSj) , on momentum and size of the electron-phonon interaction region 
in the ID case is presented in graphical form in Fig. 1. Apart from this, numerical 
calculations revealed that for each value of the electron momentum before interaction 
there is an optimum size for which the probability of phonon emission is maximal a max . 
The smaller the electron momentum, the larger is this region. 



3. Scattering amplitude 

In the preceding Section we have considered the SE for an electron interacting with 
phonons in a local nanosized region in the first approximation. Let us calculate now 
the values of the SE, E(p, e), with inclusion of the other terms of the series. We shall 
neglect, on the one hand, the effect of coherent phonon emission, and on the other, 
coherent scattering from two regions with electron-phonon interaction. This makes the 
limiting transition to consideration of a system with spatially uniform electron-phonon 
interaction impossible. It can readily be shown that E(p, e) = F(p,p,e), where the 
vertex function F(p, p', e) is given by the integral equation 

F(p, p', e) = S(p, p', e) + J E(p, fe, e)G (k, e)F(k, p\ e)dk. 

Here p,p' are the electron momenta before and after the electron-phonon interaction in 
a local region. The expression for E(p, p', e) is analogous to with the localization 
functions depending on the differences p — k and k — p'. As is well known, it is the 
vertex function that is related to the scattering amplitude / = —mFj (2tt) which plays 




Figure 1. (a) The imaginary part of the SE, plotted for the ID case vs. 
momentum p (in mc units) for different sizes of the local region expressed in units 
of do = 1/mc. In electron-phonon interaction localization regions of larger sizes, the 
dependence acquires a threshold pattern. As the size of such a region decreases, the 
threshold disappears, i.e., the probability of phonon emission is no longer zero, (b) 
The dependence of on size a for different momenta. For momenta less than mc, 
there exists an optimal size of the electron-phonon interaction localization region, for 
which the probability to emit a phonon is maximal. For momenta above mc, the size of 
the region is no longer essential. To present all the graphs on the same scale, they are 
normalized against their maximum values QT, max within a given argument interval. 



a major part in the theory of transport processes. The poles of Green's function define 
the energy spectrum of a quasi-particle, and the same applies to the / = — mF /{2ti) 
function, which is a correction to Green's function. Consider now the most interesting 
and intriguing case which is realized in the system under study. 

As already pointed out in Section 12.21 a bound "phonon" state exists in the local 
electron-phonon interaction region, but in the ID case it exists always, and in the 3D 
case - only for sizes above the critical level; therefore, the amplitude of scattering from 
state p into state p should have a pole [13J. In other words, one can write 



/( P ) = -^ (id 

K + %p 

where k, is the wave vector corresponding to the "phonon" state, and S is the correction 
to the scattering amplitude deriving from the scattering being inelastic. The magnitude 
of 6 is proportional to T,(e,p) (see ( 1101) ). For small electron momenta, 5 ~ {po) 2 - 
Significantly, it is in this feature that scattering from impurities and phonons differs 
radically from that from a local region. Indeed, as shown convincingly above, in the case 
of spatially uniform interaction 5 features a threshold behavior. For small momenta, 
the denominator of Eq. ( TTTh becomes small, with the ensuing strong growth of the 
scattering amplitude resulting from the resonance effect. For certain dimensions of the 
region of local electron-phonon interaction, the electron will undergo not only elastic 
resonant scattering from such regions but strong inelastic scattering as well. 
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4. Results 

Consider now the effects occurring only in the presence of local electron-phonon 
interaction. There are two reasons for these effects to become manifest. The first 
of them can be traced to the possibility of inelastic electron scattering from phonons 
under conditions favoring local electron-phonon interaction. For each value of electron 
momentum, there is an optimum size of the local region a max , at which the inelastic 
scattering is maximal. The second reason lies in that the creation of a bound state in a 
region with a size on the order of a c may bring about enhancement of both the elastic 
and inelastic scattering. We start with expressing the elastic, o- e iastic, an d inelastic, 
o~ineiastic, scattering cross sections through the scattering amplitude / , which can be 
done similar to the case of electron scattering from an impurity [T7] . 

^elastic = 1 — S\ (12) 
0~inelastic = ~ n(l — |5| 2 ) (13) 

where, 5 = 1 + 2ipf. For elastic scattering, we come to the relation describing resonant 
electron scattering from a bound state: 

Air 

^elastic 5 I o \ ) 

K + p z 

Examining Eq. (fT4"j) . we see that elastic scattering from a region with local electron- 
phonon interaction occurs in the same way as it does from an impurity whose potential is 
governed by the polaron binding energy and the size corresponding to the characteristic 
size of the localization region. Scattering from a local region differs substantially in 
the dependence of the inelastic scattering cross section on the electron momentum. 
As already pointed out (see Section 12. 3p . the scattering with not in this case have a 
characteristic threshold. We finally obtain that the cross section of inelastic electron 
scattering from a local region with a bound state depends on momentum as 

6(p) 2K+pS(p) 

0~inelastic — 77 9,9 \-L<-V 

p K z + p z 

In the three-dimensional case with existing bound state, this transforms for pa — > to 

^inelastic = Ana 2 - (16) 
K 

The inelastic scattering cross section can be made arbitrarily large provided the bound 
state energy is low enough, and, hence, k — > 0. 

Figure 2 plots the inelastic scattering cross section as a function of momentum 
for the one-dimensional case in which the bound state always exists. In principle, 
the dependence of a ine i ast ic on momentum reproduces on a qualitative footing the 
dependence of the imaginary part of SE on momentum. We may recall, however, that 
this system is characterized by two essentially different parameters, a max and a c . These 
sizes will become manifest in the dependence of the inelastic scattering cross section on 
size in the presence of two maxima. Although the inequality a c < a max w ao for the 
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Figure 2. Graphs of the inelastic scattering cross section ai ne i as tic for the parameter 
ratios a c /a = 1/10 ((a) and (b)) and a c /a = 10 ((c) and (d)). (a) and (c) plot the 
dependence on momentum for a given size, and (b) and (d), that on size for the specified 
momenta. As seen from (a) and (c), in the limiting cases of large sizes of the region, 
which correspond to uniform electron-phonon interaction, inelastic scattering appears 
only for momenta p > mc. For small local region sizes, the inelastic scattering cross 
section is nonzero for all, including small, values of the momentum. As evident from 
(b) and (d), for any momentum of the incident electron the inelastic scattering cross 
section features two maxima. These maxima can be identified with the characteristic 
sizes of the electron-phonon interaction region, a c and a max . The momentum p and 
the size a are presented in the same units as in Fig. [TJ To present all the graphs on the 
same scale, they are normalized against their maximum values cr ine i astic max within a 
given argument interval. 



characteristic parameters does hold, the graphs, in order to make them more revealing, 
were constructed in both cases both for a c < a max and for a c > a max . 

An important parameter in mesoscopic systems is the electron dephasing time, 
i.e., coherence loss time. At low temperatures, the dephasing time is mediated by 
various effects, among them the Kondo effect, scattering from impurities, and electron- 
electron interaction. It has recently become clear that scattering at low temperatures 
can likewise be initiated by interaction of electrons with phonons in local regions. A 
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mechanism of such interaction has been proposed for materials with clearly pronounced 
metallic properties [5]. In most semimetals, as well as in graphene and nanotubes, 
the electron gas is non-degenerate, which makes the mechanism forwarded in Ref. [5] 
inoperative. We have just shown (see Eq. (lTo*j) ) that taking into account the finite 
size of the electron-phonon interaction region also gives rise to the appearance at low 
temperatures of a contribution to the dephasing time r ph . 

In the three-dimensional case and for low momenta, the temperature dependence, 
according to Eqs. ( TT51 [TU1) . should follow the inverse pattern, r p h = 1 / npo- ine i astic ~ 1/T. 

This effect should be the strongest in one- and two-dimensional systems. The 
existence of a bound level in such systems should give rise to resonance effects, for any 
size of the electron-phonon interaction localization region (see Eq. (jlip ). Inelastic 
scattering, in particular, should appear in diamond-decorated nanotubes [TT]. An 
electron propagates through a nanotube along its axis, and regions featuring local 
electron-phonon interaction are actually regions where diamond is in contact with the 
nanotube. In this case, the dephasing time should likewise scale as the inverse of the 
inelastic scattering cross section, r ph ~ l/po"j neZastic , and the dependence of the inelastic 
scattering cross section on electron momenta and the size of decorated diamonds will 
be essentially nonmonotonic, as seen from Fig. [2] 

Graphene is a typical example of a two-dimensional system. The as-formed 
graphene can interact with the substrate. For an interaction region with the size favoring 
appearance of a bound state a c in a two-dimensional system, strong inelastic scattering 
can set in in this case as well. We believe that it is this situation that accounts for the 
observation of strong inelastic scattering reported recently [7]. 

The above analysis reveals another interesting ramification bearing on heat- 
conducting systems. We have been witnessing recently a broad research effort in the 
area of metal-insulator composites [8]. We note, however, that in a metal it is electrons 
that transfer heat, while in insulators these are phonons. Therefore, the conduction of 
heat in such composites is mediated by the effects involved in energy transfer from the 
electron to the phonon subsystem. The major contribution comes from the processes 
taking place in the surface layer. It can readily be shown that Kapitza resistance 
coefficient [H] \ at the metal-insulator interface should be inversely proportional to 
the cross section of inelastic scattering, x ~ 1 / ' ^inelastic, of electrons from phonons in 
such a near-surface region. The above consideration fully applies here. For surface- 
region thicknesses on the order of a c or a max , thermal resistance at the metal-insulator 
boundary should be substantially smaller than that in regions with other dimensions. 
This should be assigned to the strong inelastic scattering occurring in such a region (see 
Eq. fllip ). Thus, to attain effective heat transport in metal-insulator composites, one 
should provide conditions favoring creation of a near-surface region where the above 
effects would become realized. The choice of the specific size depends on temperature. 
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5. Conclusion 

The paper reports on a solution to the problem of electron interaction with phonons 
within a local region, which is based essentially on the modified Frohlich's Hamiltonian. 
It is demonstrated that such a system can support existence of a bound electron state 
within a local electron-phonon interaction region. The existence of such a bound state, 
just as that of a bound state in an impurity, translates into a substantial increase of the 
cross section of elastic scattering from such a region. 

While the probability of phonon emission by an electron in a spatially uniform case 
follows a threshold pattern, in the system under consideration the probability of phonon 
emission by an electron, even by a slow one, is nonzero. This probability of phonon 
emission by an electron will become still higher in the presence of a bound state in the 
region of local electron-phonon interaction, thus initiating strong inelastic scattering. 
Therefore, local electron-phonon interaction will affect strongly the dephasing time in 
such systems. 
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